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abstkact 


Stress distribationt have bean calonleted for a creep law predicting a 
bifher rate of plastic deformation Uian modeled in earlier studies. The expected 
reduction in stresses is obtained, altbongh quantitative results are not yet 
available because of difficulties in obtaining convergent solutions. 

Improved schemes fox calculating growth system temperature distributions 
are being evalnateil in a new subtask started at MIT. Other work in temperature 
field modeling has examined the possibility of using horizontal temperature 
gradients to influence stress distributions in ribbon. 

The defect structure of 10 cm wide ribbon grown in the cartridge system 
has been examined. A new feature is identified from an examination of 
cross-sectional micrographs. It consists of high density dislocation bands 
extending through the ribbon thickness. A four-point bending apparatus has been 
constructed for high temperature (2. 1000*0 study of the creep response of 
silicon, and will be used to generate defects for comparison with as-grown 
defects in ribbon. 

Another subtask has boen started in collaboration with the University of 
Illinois which will examine the feasibility of laser interferometric techniques 
for sheet residual stress distribution measurement. 

The swthematical formalism for calculating residual stress from changes in 
surface topology caused by an applied stress in a rectangular specimen has been 
developed, and the system for laser interferometric measurement to obtain surface 
topology data has been successfully tested on CZ silicon. 

Testing and calibration of different fiber optics temperature sensor 
configurations are underway. 


•’The JPL Flat Plate Solar Array Project is sponsored by the U.S. Department of 
Energy and forms part of the Solar Photovoltaic Conversion program to initiate a 
major effort toward the development of flat plate solar arrays. This work was 
performed for the Jet Propulsion Laboratory, California Institute of Technology 
by agreement between NASA and DOE.*' 
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I. INIEODPCTIQN 


A tatisfkctory aod^jVtliat can account for atreascs generated in eilicon 
sheet grown at high speeds is not yet available. Nuaerons atteapts to calculate 
residual stresses have been aade. but all of these suffer froa inadequacies in 
one area or another. This report describes the work in progress under this 
subcontract to atteapt to develop and test in the laboratory a stress-teaperature 
field aodel for silicon ribbon EFG. In one subtask, a coaputer code to predict 
stress-*teaperatnre field relationships in steady-state sheet growth is being 
developed at Harvard University. The stress state is paraaeterized by a 
two-diaensional teaperature field and growth speed. Incorporstion of time 
dependent stress relaxation effects is through a creep law to aodel the iapa'ct of 
plastic flow on the sheet residual stress state. A second aspect of the program 
deals with the development of a aodel to predict the teaperature field in a 
moving sheet from given system component temperatures (i.e., the sheet 
environaent) , and studies ezperiaental aeans to verify the model. 

The results froa modeling of stress in silicon sheet with creep show that 
plastic deforaation processes can significantly reduce stress during steady-state 
high speed sheet growth. The initial atteapts to compare calculated stress 
distributions and experiaental results in the case of growth in the 10 ca EFG 
cartridge system raise a number of iaportant issues. These relate to the 
directioLS for future research needed to provide a test for the aodel and 
inforaation relevant for achieving reduced stress growth configurations in 
practice. These issues are examined here. 
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The ialtiel eneljtis hai invest Ageted the iapact of two levels of creep 
intensity. The creep response hfis been approxiasted by choosing a single 
exponent (n “ 5) for the stress dependence of the strain rate and a teaperature 
independent activation energy (~S eV) for defect action. This has its 

liaitations. as evidenced by data froa varions workers, which show that these two 

paraaeters are a strong function of stress aagnitnde and of teaperatnre. In a 
second approxiaation, only steady-*state creep has been aodeled, and transient 

neglected, lik view of these approxiaations, it is not possible to choose 
a repr;ns«iitative experiaental datna point for a fit of the constitutive law used, 
as has been done up until now, and expect to arrive at an understanding of creep 
influence and at a test of the aodel at a quantitative level. Additional data on 
creep response relevant to silicon sheet growth are required to gain 

understanding both of the effect of iapurities, such as carbon, on defect 
generation and notion in ribbon at high temperatures, and of relaxation phenomena 
at the stress and strain rate levels estimated on the basis of the present 
calculations. A program on study of high temperature creep has been started to 
provide more information in these areas. 

The sheet temperature profiles u od in the stress analysis to date all 
lead to large elastic stresses, sufficient to result in buckling and ribbon 
fracture. This includes the profile calculated for ribbon grown in the 10 cm EFG 
cartridge system. The residual (room temperature) stress calculated with creep 
included is sensitive to the details of the temperature profile in the active 
afterheater region of the cartridge. However, the stresses within about 5 mm of 
the growth interface, associated here with high temperature buckling, do not vary 
significantly and are not reduced appreciably for the modifications in the 
cartridge temperature profile examined. Attempts to vary growth conditions in 
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the oertrldge ia order to prodnce ribbon with different ttreei levels to provide 
a test of the nodel have net with linited snccess. They show that the 10 on 
cartridfe generally laehs the flexibility necessary for this work. Consideration 
is being given to design of a ribbon growth aysteai more suitable for this 
purpose . 

Calculated sheet temperature profiles have not yet been confirmed by 
measurements in the ribbon. Development of a temperature sensor based on fiber 
optics, which has the potential to provide c^ata at a neceasazy level of accuracy, 
is in progress. This confirmation of calcnla;ted profiles is desirable because 
details of the heat transfer in the crucial region near the growth interface need 
to be ascertained. The modeling shows that the temperature profile shape is very 
sensitive to values used for several silicon material parameters that are not 
well known. This adds an uncertainty to the validity of the calculations beyond 
what arises from approximations used in the heat transport model. Refinement of 
the temperature profile calculations is in progress through application of more 
advanced finite element schemes for modeling of heat transport in the meniscus 
melt and ribbon. 

Average residual stress levels have been estimated in 10 cm wide ribbon 
grown with several different temperature profiles using a scribe and split 
technique. These are found to be in the range of 10 to 30 MPa, and generally are 
an order of magnitude or more below the maximum predicted sheet elastic stress. 
The modeling shows these low measured stress levels will only be achieved through 
creep. The technique used cannot provide information on the stress distribution 
in the ribbon. Several other approaches for measuring strc are being 
evaluated. These include multi-finger cutting and the use of laser 
interferometry. An improvement in the sensitivity of the stress measurement 
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tdol' tits* !• partieultrly aeoetsary la order to •valoate the iapaot of edge 
effects oa the orersll stress distribatioa. These arise froa thiokaess 
aoaaalforaity at the ribboa edge, created priaarily by growth froa balboas dies 
as ad to eahaace edge stability. Deviatioas la both residaal stress distribatioa 
aad high teaperatare creep aad backliag froa those predicted for a aaifora 
thiokaess ribboa caa be expected to arise. 

The goal of the prograa is to coabiae the resalts of the above areas of 
stady to arrive at a aodel that caa predict stress-teaperatare field 
relatioaships ia steady-state siliooa sheet growth aader realistic ooaditioas. 
Miaiaaa stress ooaf igaratioas will be soaght, aad aa atteapt to ooastract aa EFG 
siliooa ribbon growth systea that caa verify this aodei will be aade if it 
appears saoh ooaf igaratioas caa be achieved ezperiaeatally. 


II. PKOGBESS KEPOtT 

Th* •xperiatntil work on this progrui hat been reetrnotnred into four 
subttski (see also the revised Prograa Plan in Appendix I): (A) Advanced Systen 

Design. (B) Creep Law Studies, (C) Optical Fiber Te^eratnre Sensor Studies, and 
(D) Eeaidnal Stress Neasnrenents. The theoretical work on stress analysis will 
continue to study reduced stress configurations in coordination with subtask (A) 
above . 

fork in progress in these areas is described in detail below. 

A. StiUb fisiAia 

1. Stress Analysis (J.f. Hutchinson and J.C. Laabropoulos, Harvard 
University) 

Calculations of stress distributions have been carried out to 
investigate the effect of changes in the constitutive relation for creep from 
that used in earlier studies. This is obtained froa new experiaental data,^ 
which give a strain rate dependence on stress and temperature considerably 
different froa that used in the stress analysis up to this tiaa. This relation 
is obtained froa work hardening data taken between lOOO'C and 1300*C, and 
encompasses higher (*>10-* s-*) strain rates than used in earlier studies. 

The new creep law has the fora: 

|C . Cf(T) 

with C - 5.85 X 10»» (OPa-s)-*, f(T) - [exp(-41800/T)/T] , and n - 3.6. 

CalcuJations have been carried out for a speed of 6 ca/ain for the cartridge 
system temperature profile. The results of the calculations are shown in Figs. 1 
to 3. Complete convergence of the solutions was not obtained, and so the results 
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Fig. 1. Comparison of variations across ribbon width at room temperature (x = 20 cm) 
at V = 6 cm/min for two creep laws: 1. high creep condition, 2. new creep law. 
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Fig. 2. Comparison of Oyy along ribbon centerline at V = 6 cm/min for; 
1. high creep condition, 2. new creep law. 
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line at V = 6 cm/mln for: 




oAn oaljr b« viAV«d oa a qaalitativa laYal. Thla aat dua to aa iaadaqaaoy in tba 
HXid apaolat ia tha raeloa of blfb atraia rataa aaar tba iatarfaaa. Tba 
oalealatioa grid will bava to ba ravorkad la ordar to alialaata tbia problaa. 

Kadmatioaa ia o^ aad both ooear, aa aapa<jtad, vitb tbia aav 

eoaatitatlTa ralatier., wbicb pradiota bigbar atraia rataa (Fig. 9) tbaa foaad for 

araa tba "bigb" eraap eoaditioa atadiaa la oar iaitial aaaljaia. o baa ita 

TT 

paak Talaaa radacad aboat a factor of two. aad tba rooa taaparatara atraaa walaaa 

alao ara radacad by aboat a factor of two. Tba dlatribatloa of tba lattar aoroaa 

tba rabboa width baa a caatral ragioa of taaaila atraaaaa wbicb gradaally 

daoraaaaa aad baeoaaa coaq^rcaaiwa at tba ribboa adga (Fii|. 2). At lavala of 10 
MPa. tba rooB taaparatara atraaaaa ara comparable vo tboaa aatiaatad by applyiag 
tba acriba aad aplit taebaiqaa to 10 ca wida ribboa grown with tba ca*-tridga 
ayataa aodalad. Tba paak atraia rataa ara factora of two graatar (Fig. 3) and 
aigaificaat ragioaa of aoa-*xaro atraia axtaad to tba lower taaparr^mra ragioa of 
tba aftarbaatar. 

2. filAi A. Brown aad E.M. Bttonaay. MIT) 

Fiaita alaaant aaalyaia of tba aiaaltaaeona affacta of beat 
traaafar and capillarity baa extaada> ^raviona work* to inoorporata aa affieicat 
Nawtoa iteration acbaaa aad to iaolada tba dataila of tba fluid flow in tba welt. 
Ika Newton acbaaa lead to a factor of tbraa aawinga in ealenlatioa *’.iaa and to 
tba iaplaaantatioa jf eoapntar-aidad aatboda for tracking aolatioa faailiac ia 
taraa of a (iron paraaatar. a.g.» growth walocity aat-poiat taaparatara or 
atatic head b^^^. New operating liaita for incraaaiag tba atatic bead on 
predicted. Tbaaa calcalatioaa ara wary algnificant. Wban taken togatbar with 
tboaa reported in Bttonaay at al.* tba finite alaaaat analyaia daacribaa tbo 
"operating region" in (V^. T^, apaoa for a particular die gaoaatry aad 
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tmhitnt heat transfer eaviroBaent. Ezprestinc the operating liaita as V • 

f 

gives the aaziana aecessible growth velocity for any particular 
coabination of and £,££• Ve are exploring nov the effect of changes in 
capillary width on this growth rate ceiling. 

A direct oMaparison has been eoapleted between calculations which include 
the details of the velocity field in the aelt and ones using a Modified unifora 
velocity field that accounts for the tapering of the aeniscus. The two 
calculations were found to be in good agreeaent in silicon for growth speeds 
below S ca/ain. The detiils of convective heat transfer becoae iaportant in 
silicon only near growth speeds of SO ca/ain. but at auch lower speeds for the 
solidification of aaterials with lower themal conductivity, e.g.. Al^O,. 

A coaparison betw^vu aeasureaents of sheet thickness as a function of 

growth speed is being carried out. Using the aabient teaperature distribution 

aeasured by Ealejs and Bell it Mobil Solar, thickness versus growth rate data* 

will be fit over the entire range of pull speed using the die set-point and 

die aabient teaperatures as adjustable parameters. Of these, the die 

set-point value is the aost iaportant. Plans have been aade for a systematic set 

of ezperiaents with varying ^ ^ and V for a acre coaprehensive check of the 

0 eti g 

EFG Model and the operating liaits predicted by it. 

3. EFG Svstea Teaperature Field Modelina (R.O. Bell and J.P. Ealejs. 

Mobil Solar) 

It has been suggested that horizontal teaperature gradients 
(across the ribbon width) nay be used to coapensate for the effects of high 
vertical gradients and to reduce the overall level of stress present during sheet 
growth. We have investigated the possibility of doing this in the 10 ca EFG 
systea by introducing horizontal heat flux variations into the teaperature field 
aodel and calculating resulting horizontal ribbon isotheras. 
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I& Bost of tho toaporatnre aodeliag oolculatioiit, the ribbon tnrfaoe heat 

flux haa been calculated ascuaing that the horixontal teaperatuxe was constant 

across the width of the ribbon enclosure, beyond which (i.e., at the enclosure 

end wall) it was rooa teaperatnre (300*K) . Thus the flux is essentially constant 

until very near the ribbon edge. To siaulate the effect of an enclosure wall 

horizontal gradient, the ribbon surface flux was assuaed to decrease as if its 

effectiwe teaperature decreased 100*C froa the center to the edge with a 

quadratic ourwature, i.e., T ■ T . - 100 (y/W)*. 

orx coiibor 

Figure 4 shows oaloulated horizontal ribbon teaperatures snperiaposed on 
the profile along the center of the ribbon. In all cases the horizontal gradient 
is zero at the ribbon edge because of the boundary conditions. The thoraal 
conductivity of the silicon saooths out the horizontal variations so that the 
largest edge to center teaperature difference in the ribbon is 55*C rather than 
the 100*C iaposed. The aaziaua horizontal gradient is of the order of 35'C/ca 
for this 300 (la thick ribbon. This is very saall coapared to the 1500*C/ca 
vertical gradients nea the interface but coaparable to those 3 or 4 ca froa the 
interface. 

The large interface gradients will aake it difficult to tailor the 
horizontal teaperature to coapensate for stresses produced near the interface, 
but it aight be possible to deal with stresses generated a centiaeter or so away 
froa the interface. The question still reaains as to which region doainates the 
rooa teaperatnre stress distribution. This will be strongly dependent ca the 
creep rate and its teaperature variation and on the detailed shape of the 
teaperature profile. Creep experiaents are in progress at Mobil Solar and are 
discussed further below. 
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Teaperituro field deti were inoospletely reported in the third quarterly 
report^ became of the oaiaaion of one oaae froa a table. Thus. Table II on page 
16 of thia report it reproduced here with a fourth cate. Cate D. included. 

Table II. Theraal Conditiona and Geoaetry uted in Modeling. 

(Corrected Table froa Third Quarterly Report.) 



(hirve A 

B 

C 

D 

Ribbon Thickneaa 

300 pa 

300 pa 

300 pa 

300 pa 

Bate Plate Teaperature 

1500*C 

1200* C 

1200«C 

1500»C 

Cold Shoe Teaperature 

4S0 

450 

400 

400 

Afterheater Teaperature 

1200 

1200 

960 

960 


1310 

1370 

1490 

1427 

On account of this oaiaaion, wrong 

identification of the 

temperature 

diatributiona used in atreaa 

calculationa 

war Bade. 

Curve A 

ahould now 

correapond to the EFG profile. 

Caae 1, and 

Curve D to 

the reduced 

afterheater 


profile, Cate 3, of Fig. 1 of the third quarterly report. 

B. Creep Law Studiet 

1. Four-Point Bendina 

A four-point bending apparatue hat been conatructed for uae at 
high teaperaturea to atudy creep in ailicon and it at preaent being teated. It 
it aade of graphite in order to withatand teaperaturea up to the aelting point of 
ailicon, and followa fairly cloaely the ASTM E328 recoauaended apparatua. A aide 
view of the apparatua it ahown in Fig. 5. 

The firat ezperiaenta have been carried out on initially defect-free CZ 
ailicon wafera and on Silao ailicon. Saaplea were out into approziaately 5 oa by 
10 oa blanka. CZ and Silao aaaple thickneaaea were 0.3 and 0.4 aai. 


13 



OKIGINAL PAGE 13 
OF POOR QUALITY 



Fig. 5. Four-point bending apparatus for 
high temperature experiments. 



Bxporlaants have bean oarriad out at two tanparaturaa, 121S*C and 1360*C. 
Tba taaplaa were atraaaad nndar prograssivaly inoraaaing loads nntil a total 
daflvotion of about 0.10 to 0.12 oa was obtained. Than the load was raaoved and 
the saaplas w9T» cooled to roosi taaparatnxe. The displaoesiant vs. tiaa graphs 
are given in Fig. 6. 

The idealized stress distributions along the saaple length and through its 
thickness are shown in Fig. 7 for the four-point bending configuration used here. 
The aaziaua surface stress (tension on the bottoa and ooapression on the top) is 
calculated froa the aoaent of inertia I and bending aoaent M as: 

NC 

«s ■ r 

I - 7T ft* M - FI 

12 

For our ezperiaents W * 5 ca and I ^ 1.9 ca. Thus for loads of 100 to 200g, the 
aaziaua values are of the order of 10 MPa for typical saaple diaensions. For 
the data of Fig. 6, values up to about 15 MPa were realized. 

Future plans for the creep ezperiaents are to study the defect structure 
generated in CZ wafers as a function of loading conditions, and to ooapare it to 
as-grown ribbon structure. 

2. Ribbon Defect Studies 

Ribbon cross sections have been ezaained to try to identify defect 
features and their distribution which aay give inforaation on stresses acting 
during ribbon growth. The ribbon was sectioned perpendicular to the growth 
direction, polished and etched to reveal the cross-sectional structure. Several 
tjrpical aicrographs are shown in Fig. 8 and 9. A new class of dislocation 
arrangeaents is shown there. These occur in the fora of narrow bands of very 
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Fig. 6. Displacement vs. time graphs for (a) CZ, and (b) Silso 
silicon wafers stressed by four-point bending. 
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Fig, 7. Geometric and stress relationships for 
sample of thickness t (and c = t/2) in 
four-point bending. 
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. Cross-sectional micrographs of 10 cm wide ribbon grown at 2,6 cm/min 
and with afterheater setting of 1000°C. (a) Magnification X70.4, 

ribbon thickness is 0.55 mm; (b) magnification X45.2, thickness 0.15 m. 
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Cross-sectional micrographs of 10 cm wide ribbon grown at 
3.5 cm/mi n and afterheater setting of 1100 C. (a) Ma9- 
nification is X70.4, ribbon thickness is 0.28 mm; (b) high 
magnification (X652) view of area identified by arrow 
I (a). 





high disloottion density whloh extend thronih the orost section. These ere 
typioelly of the order of the ribbon thickness in extent, end ere bounded by whet 
eppeer to be slip bends or yet sore intense dislocetion bends. Fignre 8(e) shows 
whet eppeers to be en eerly stege of the bend dewelopeient, where the dislocetion 
density is intemediete between everege levels occurring in edj scent regions end 
thet in Bore highly stressed bends. Exeaples of the letter ere found in Fig. 
8(b) . The bends occur often in peirs or triplets sepereted by distences of the 
order of severel tines the ribbon thickness. Double bends, such es shown on the 
right of Fig. 8(b) ere el so ooBBonly observed. 

A high Bsgnif icetion aicrogreph of one of the more intensely disloceted 
bend edges is shown in Fig. 9(b). These regions do not eppeer to be oriented 
elong the twin plene directions (one of which l.^.es elso in this aicrogreph), but 
rether eppeer to aeke lerger engles with respect to the surfece then do the 
typicel twin plenes seen in the cross sections. 

c. £i^ QsHsm. TfWfUtyrt gfMPr 

A probe wes constructed to hold the high end the low teaperetnre 
portions of the fiber optics sensor in e oonf iguretion thet would ellow it to be 
pieced conveniently into the efterhester region of the 10 oa certridge. The 
first test of the probe in the certridge wes to aeesnre teapereture in the 
hottest pert of the efterhester, end to use the efterhester theraooouple es e 
reference. Initiel tests of the probe showed the signel le\el wes aueh higher 
then thet occurring in the test furnece used for the initiel sensitivity studies. 
It wes decided to go to e two wevelength aeesnreaent in order to aeke ebsolnte 
teapereture deterainetions end so evoid the probleas essooieted with verying 
light intensity levels. 
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Tht two war«ltii|tlit ehosaa for tlio MosaroaoBt woro O.tf tad 0.7 |ui. Tho 
probo waa oalibratad at roM toaparatara with a apaotral radioaatar aaiag a 
toafataa ELB laap oparatiag at a taaparatara of 3320*K. Tba aftarbaatar ragioa 
taaparataraa wara apaia aoaitorad. It waa foaad that tha proba ealibratioa 
factor obtaiaad at rooa taaparatara atill did aot fira taaparataraa that 
eorraapoadad to tba aftarbaatar tbaraocoapla raadiaga. 

It appaara that tba aajor problaa with tba proba at praaaat ia that too 
aacb atray radiatioa ia baiaf aceaptad iato tba fibar optica circuit that 
orifiaatac froa ragioaa otbar tbaa tba proba tip. A aacoad probe ia baiag 
ooaatruotad that will ba placad ia a cooler ragioa of tba cartridge aaar tba cold 
aboaa to attaaipt to raduca tba affaota of atray radiatioa. 

0. Raaidual Straaa Maaaaraaaata (A.T. Aadoaiaa aad S. Daaylnk, Uaivartity 
of Illiaoia at Chicago) 

It baa racaatly baaa propoaad that laaar iatarf eroaatry way ba oaad to 
aTaluata raaidual atraaa diatributioaa on a aacroaeopic aeala la abaata of a 
■atarial aueb aa ailieoa. la tbia taebaiqua. tba laaar light la eaad to gaaarata 
aurfaoa topograpba of tba abaat aurfaca uadar diffaraat atataa of applied atraaa, 
aad tba daforaatioaa aaaaurad are uaad to back calculate tba raaidual atraaa ia 
tba abaat. Tba aatbauatical foraatioa for calculatiag tba raaidual atraaaaa baa 
baaa davalopad, aad ia giraa for a raotaagular aaapla gaoaatry ia Appaadix II. 
Tba apparatua for tba topology aaaauraaaat baa baaa coaatrnotad aad taatad to 
aucoaaafully back calculate aa applied ia-plaaa atraaa ia a CZ ailieoa wafer. 
Thaaa raaulta are diacuaaad ia Appaadix III. Tba aaxt atap ia to attaapt to 
■aaaura raaidual atraaa ia EFG ribboa. 
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ABSTRACT 


Residual stresses in silicon ribbon can result from a two-dimensional 
temperature field or a change in the growth speed of the ribbon and the 
magnitude and sign across the width of the ribbon can vary. These stresses 
can be related analytically to an intentionally applied out-of-plane bending 
moment and the curvature of the distorted ribbon, quantities which can be 
easily obtained by a unique experimental measurement. This report is a 
summary of an analysis assuming two conditions on residual stresses in 
silicon ribbon: in-plane tension and side shear and in-plane compression 
and side shear. 

INTRODDCTION 

Residual stresses are extremely important in silicon ribbon since these 
stresses can cause buckling and in some cases fracture of the ribbon. A 
knowledge of the magnitude and sign of the residual stresses can therefore 
be used to optimize the growth conditions of silicon ribbon. However, at 
present, our knowledge of residual stresss in silicon ribbon is sketchy. 
Some recent modeling has shown [10] that thermal gradients and change in 
growth speed of the ribbon can result in non-uniform residual stresses 
across the width of the ribbon. Experiments are underway to verify these 
stresses but these experiments are destructive. 

Several non-destructive techniques are available for determining 
residual stresses. The two most widely used are the X-ray technique and the 
hole-drilling strain gauge method [1-7]. Both techniques have limitations. 
The X-ray technique is limited since X-ray absorption occurs within several 
micrometera of the surface depth. The resolution of this technique is 


2 

'\i 5 kti over teverel ma* in eree. In addition, stresses throughout the 
thickness eannot he detemined conveniently. Another technique is the hole- 
drilling strain gauge siethod where a special rosette strain gauge is bonded 
to the ribbon and connected to a strain aeasuring instrunent. A small hole 
is then drilled through the center of the gauge and the resulting change in 
strain around the hole due to relaxation is measured. This technique is 
semi-destructive and measures the strains within the top surface layers 
(usually a few millimeters) and a very specific location. A determination 
of stresses is obtained from the strains. Both techniques are 
unsatisfactory for measuring residual stresses in ribbon because (1) the 
thickness of silicon ribbon is '^<1 mm or less and stresses through the 
thickness are desired ?n. (2) residual stresses over a large area are 
needed. 

Residual stresses can be determined by double beam interferometry which 
relies on the application of flexural loads and the measurement of the 
deflection and curvature of the silicon ribbon. Out-of-plane deflection 
will either be accentuated or reduced depending on the sign and magnitude of 
the residual stresses. The residual stresses are then determined by 
comparing the measured and predicted deflections and curvature. Any 
difference is ascribed to residual stress. 

This technique has a number of advantages over existing non-destructive 
techniques for the measurement of residual stresses f.l~7] . First, large 
spatial areas can be examined in contrast to the X-ray or blind hole 
drilling techniques. Second, through-the-thickness measurements are made as 
opposed to only the top surface layers as in the other mentioned techniques. 
One disadvantage is that the residual stresses may vary significantly 
through-the-thickness and the interferometry technique will average these 


values . 


In this progress report, we present the snelysis for the relstionship 
between the deflection and curvature and the bending moment and residual 
stress in thin rectangular ribbon. This analysis will be used in 
conjunction with experimental measurements on determining the residual 
stresses . 

ANALYSIS 

We make use of a stress analysis by Hutchinson and Lambropoulos [lOj of 
steady-state ribbon growth. Their model takes into account a two- 
dimensional temperature field and steady state growth conditions of the 
ribbon. After cooling, the stress conditions are 

-30 MPA < a - f(y) < 30 MPa 


where x is in the direction of ribbon growth and y is along the width. 

Figure 1 shows typical residual stress variation along the width 
of the silicon ribbon tepresentatlve of those found [lO]. Due to the 
heterogeneous nature of a, residual shear stresses will be produced in the 
x-y plane. An analysis of the residual stresses in the ribbon must include 
the tensile, compressive and shear stresses in the ribbon. Figure 2 shows 
qualitatively the variation of shear flow, q, along a longitudinal section 
(parallel to x-y plane). This variation in 'q' can be assumed to be linear 
except at the end sections where the shear stress should be zero since the 
ends of the ribbon are traction free. The assumed variation in shear flow 
shown in Fig. 3 can be expressed as 

( 1 ) 

where E represents the change in shear flow per unit length. 
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The verietion of along the width identified in Fig. 2 necessitates 
the nae of separate models for different states of stress. Formulations for 
two different oasea are presented in the following sequel. 

Case I: Rectanaular Plate Under Bendina . In-Plane Tension and Side Shear 

If the ribbon shown in Fig. 2 was subjected to an out-of-plane bending, 
the mid-portion could be treated as a rectangular plate under bending, in- 
plane tension and side shear. The state of loading for a narrow strip of 
width is shown in Fig. 4. 

If the deflection along the ribbon is approximately as 

w - Ax^ + Bx + C 


and the boundary conditions are given as 
[x “ 0, w = 0] 

[x B L, w “ 01 

then 

w - A(x^ - Lx) (2) 

Considering angular equilibrium. 


where 


M- - M - P. . M, 


r 


qw dx 


(3) 


(4) 


From linear equilibrium the stress residual in the x— direction can be 
written as 


P - 2 



q dx 


( 5 ) 


substituting Eq. (1) into Eq. (5) 


P - 2 


ORIGINAL PAGE IS 
OF POOR QUALITY 


solving for E 

P 



and comparing Eqs. (1) and (6) we get 



Sobstitnting Eqs. (9) and (7) into Eq. (11) 



Snbstitnting Eq. (8) into Eq. (3) 

H- . M - P, . - Lx„) 


( 6 ) 


( 7 ) 


(8) 


(9) 


Finally assoming small deflections and utilizing the relation for plate 
cnrvatnre, we get 

1 _ d^w ^ _ M' 

p “ dj? “ D 

where 

^ Eh^b 
” 12(1-V^) 

h • plate thickness 
b ■ plate width 
E > modnlns of elasticity 
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V ■ Poisson’s ratio 

The governing equation for Case I can then be written as 






where 



•2k (to be detemined experimentally) 


N “ applied moment (known) 
w “ measnred deflection at x^ 


( 10 ) 


Next a portion of the ribbon subjected to compressive residual stresses 
(between A and C or G and I) will be considered. 


Case II: Rectanaular Plate Onder Bending, In-Plane Compression and 

Side Shear 


The state of loading for this particular case is shown in Fig. 5. 
Similar to Case I, the angular equilibrium can be expressed as 


M' - M + Pw - 


( 11 ) 


Substituting Eq. (8) into Eq. (11) 

M- - M . P, - ^ (xj - U„) 


and the governing equation 


d^w HP AP / 2 

- 5 ■ 5 ’ ^ r (*o 



( 12 ) 


(13) 


As can be seen from Eqs. (10) and (13), the curvature at any length 
along a deflected ribbon is dependent on M, the applied moment and P, the 


stresa residual. 


To determine the dietribution of residual stresses a moment is applied 
and the deflection is recorded. A parabolic curve fitting procedure is then 
used to fair parabolas through a series of data points along z and 

j axes in order to find the slope and curvature of the deflected surface 
at desired locations. Finally, experimentally-computed values of slope and 
deflection can be need in the theoretical model presented here to find the 
in-plane residual stress distribution. 



y 
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Fig. 1. Schematic representation of residual stresses along an arbitrary 
transverse section of a silicon ribbon. 
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Fig. 2, A schematic representation of the res‘'dual stresses along the length 
of a silicon ribbon and the shear fli on longitudinal sections where 
the sign of residual stress changes. 



Fig. 3 Shear flow variation along a longitudinal cut in a silicon ribbon. 



M « (applied moment per unit length). Ay 

M'= internal bending moment at x . 

‘ (stress residual along x per° 
unit width). Ay 


Fig. 4 State of loading in Case I 



Fig. 5 State of stress in Case II. 

( M, M' , P are as described in Fig. 4) 
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ABSTRACT 

At later interferoaetry technique it need to determine 
dcflectiont of tingle oryttal theet tilioon to an accuracy of 316 na. 
Theae deflectiont are ntilized with an analytit to yield the til icon 
through-thicknett retidoal ttretiet. The tentitivity of the technique 
to extract reeldna atrettet haa been determined by application of in- 
plane loadt that timnlate the retidnal atrettet and uting the plate 
deflection and onr theoretical analytit to calculate the in-plane 
loadt. The calonlated and applied in-plane loadt corretpouded to 
within 8% for atrettet of the order of 100 kPa (14 pti). 


i wm ppcTi w 


Besldvil ttratttt ar« «xtrea«ly laportant in ailicon ribbon since 
these stresses can oanee bnekling and in some casea fracture. The 

importance of the.# streaaea has been ellaoidated in pre\Ious reports, and 
a number of workers are attempting to calculate and measure these residual 
stresses as a function of growth parameters of the ribbon. We are 
experimenting with a laser interferometry (non-deet motive) technique that 
utilixes a conventional laser interferometry apparatus In conjunction with 
a novel loading geometry and analysis to extract residual stresses in thin 
silicon sheet. Since the technique is new. we have attempted to determine 
its sensitivity by testing rectangular-shaped single crystal Cxoehralski 
silicon whose residual streases are assumed to be zero. We have applied 
in-plane loads which siamlate residual stresses and we have used the 
experimental data and analysis to extract these loads. 

EIPERIIIENTS 

A block diagram of the essential elements of the technique and a 
schematic of the experinent are shown in Figs. 1 and 2. A single crystal - 
0.6 am thick Cxoehralski silicon wafer as scribed to have a rectangular 
shape. The specimen was flexed under 4-point bendioi to produce a constr.nt 
bending moment in the mid-portion of the plate. A known i»-plane load P 
that simulates residual stresses was applied and the resulting interference 
pattern under different load levels (tensile and compressive) was recorded 
photographically. 

The mid-point deflection and plate curvature was found to increase 
under compressive loads and decrease under tensile loads. This situation 
it depicted schematically in Fig. 3. A schematic of the top view of the 
silicon sample is given in Pig. 4 where the inner 


load bearing pins of the 4-point bend apparatus (2-in span), the applied 
load P and noment the line along which the optical data were taken and 
the coordinate frame (x and y axes) are shown. The local deflection and 
cnrvature was found by analyzing the interference pattern and by using 
geometric scale factors (to determine the actual distances from the 
photographs). The interference pattern for five combinations of P and M 
are shown in Fig. 5. As can be visually identified, the relative positions 
and density of fringes vary with P and N. 

ANALYSIS 

(1) Determination of the Local Curvature From the Experimental Data 

To find the local curvature of the silicon plate at an arbitrary point 
Q, 3-point method is used in conjunction with the optical data in the 
vicinity of point Q, Fig. 6 shows a small portion of the silicon plate 
which contains point Q at which optical information is gathered. 

The radius of curvature of the silicon plate at point Q is obtained 
by passing a circle through! the data points which has the form 
+ 2dx + 2sZ + f = 0 (1) 

and the radius of which can be expressed as 

R =(d^ + e^ + f ) (2) 

where 

d, e and f are constants to be determined. 

After peducicg the optical information at three solution stations 
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surrounding point Q we get 

Xj Zj Kt point 1 , 

X 2 Z 2 at point 2 and 

Xj at point 3 

Su.bstituting these data into Eq. (1) we get. 

2 dxj + 2 eZ 2 + f » - + Z^) 

2 dxj + ZeZj + f = - (xj^ + Z 2 ^) 

2 dx 3 + ZeZj + f - - (xj^ + Zj^) 

or in matrix form 


xi Zi 

*2 
X 3 Z 


Z 2 1 


2 d 

2 e 

f 


xj- + Zj 


*2 * ^2 


*3 * h' 


from which d, e and f can be evaluated to determine the local curvature R. 
(2) Determination of for a Given Bending Moment M and Axial Force P . 

Fig. 7 shows a portion of the rectangular silicon plate under the 
applied loading. Assuming the maximum deflection to be less than the plate 
thickness, the radius of curvature of the convex surface at point Q can be 
written as 


h Ph^ + Ebh^ 

•*o = T — nr — 


where h = plate thickness, 

b = width, 

E = Young's Modulus, 
M' = M - Pwq 


P = in plane axial load. 


tf = applied bending moment and 
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Wq = deflection of silicon plate at point Q. 

The radios of curvature at the midspan is then given as 

, h Ebh^ Ph^ 

•>o|, -I 

RESULTS 

The sensitivity of the technique is obtained by measuring 
obtaining from the optical data, calculating P and finding the 

difference between the applied load and the experimentally computed load. 
The smallest in-plane load used was 1.982N which results in a stress of 
77.24 kPa ( 10 psi) . The results of the percent difference between the 
calculated and computed in-plane lead for different loading conditions are 
given in Table I. The percent difference ranges from 0.036 to 8.26 
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Figure 1. 





ORIGINAL PAGE IS 
OF POOR QUALITY 




incident reflected 
Dork Fringes 
Destructive Interference 


if 2j0 = n X 



incident reflected 
Light Fringes 
Constructive Interferenn 


Figure 2 
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Figure 4. 





0.02685 N.m 


P = 3.934 N 


M = 0.0403 N.m P = 3.924 N 
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u.uituj N.m = -i.yoz in 


Figure 5. 
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Figure 7. 
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R (m) 

cSlculated frm 
optical data 

M (N-m^) 
applied 

P (N) 
applied 

P (N) 

calculated 

Difference 

1 

mm 

6.392 _ 


BHilli 


HiSBIIIi 


6.3914 

0.02685 

1.962 

1 .8 

8.26 

mmmt 

6.3938 

0.02685 


3.8 

3.2 

mmmm 

4.2903 

ISRIBSSH 


3.82 ... 



m 

4.243 

0.0405 

•1.962 

-1.88 

m 


* (percent difference calculated in bending moment M ) 


Table I. Percent difference between calciilnted and applied 
in*plane load , P , for a given bending moment M. 



















